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*  Goals 

*  Review  of  Past  Work 

*  Argon  Collisional-Radiative;  comparison  of  reduction  mechanisms 

*  Non-Maxwellian  CR 
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Goals 


Utilize  hybridization  techniques  to  produce  accurate  and  efficient  plasma 
simulations  that  spans  many  orders  of  magnitude  in  both  space  and  time. 

Capture  complex  physics:  excitation/ionization,  transport,  radiation,  etc. 

Consistent  collision  operator  across  different  levels  of  fidelity. 


Current  Focus: 

•  Generalization  of  collisional- 
radiative  kinetics  with  level 
grouping 

•  General  Hybridization  techniques 

•  Focus  on  each  solver  before 
hybridization 

•  Special  attention  to  low  density 
low  energy  conditions 


Laser  Plasma  Interaction 
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Summary  of  Past  Work 


Maxwellian  Inelastic  Collisions 


•  Detailed  CR  model  for  multiple  ionization  stages 

•  Validation  against  experimental  data 

•  Nonequilibrium  radiation  transport:  coupling  with  a 
collisional-radiative  model 

•  Inelastic  collisions  in  a  MF  plasma:  enhanced 
thermochemical  kinetics. 


Multiscale  Hybridization 


•  A  time-parallel/multiscale  method  with  energy 
preservation 


Boltzmann  -  Uniform  Grouping  comparison 


•  Boltzmann  is  more  accurate  when  compared  to  electron 
configuration 


Computation  Time  (sec) 
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Collisional  Radiative  (CR)  Overview 


Updates 


•  Investigated  Quasi-Steady-State 

•  Investigated  single  temperature 
Boltzmann  grouping 

•  Reintroduced  radiation  sources 

•  Investigate  grouping  sensitivity 

•  Linked  with  LANL  database  for  Argon 
cross  sections  and  atomic  level 
information 

•  Algorithms  not  hard  coded  for  Argon. 

•  Adaptive  time  steps 


Kinetics 


•  Electron-Impact  Rates: 

•  Excitation  <->  De-excitation 

•  Ionization  <->  Three-body 
Recombination 

•  Radiation  (Planckian) 

•  Photoexcitation  <->  Stimulated  + 
Spontaneous  Emission 

•  Photoionization  <->  Simulated 
Radiative  Recombination 
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CR  Grouping  Techniques 


Quasi  Stead-State  (QSS) 


•  Assumes  fast  kinetics  between 
states  within  an  ion  distribution 

•  Assumes  longer  diffusion/decay 
times  than  excited  state  lifetimes 

•  Does  not  conserve  mass  of 
excited  states 

DN?k  ^  _ 

±  D  +  fc  I7  +  k  i\T+k 

~  n  *n  fVn 

DNfk 

— — —  =  0  for  n  >  1 
Dt  ’ 


Boltzmann  Grouping 


•  Conserves  number  density 

•  Preserves  energy  in  groups  through 
group  temperature  description 


•  Consistent  with  Boltzmann 
equilibrium. 


Conserved  Variable: 

Nr 


Nno&K  = 


n0 


9n0  j 


I 


9iei 


^Ei/Tn 
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Effective  rate: 


g.e  AEi/Tn 


i 
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Complexity  Reduction  for  Argon 


Prevalent  in  numerous 
applications,  i.e.  electric 
propulsion,  LPI,  etc. 
Extracted  atomic  data  sets 
from  Los  Alamos  Natl.  Lab 
(LANL)  database  for  all  18 
iso-sequences  of  argon 
Detail  level  of  atomic  data 
is  LS-coupling.  Reduced  the 
data  to  detailed 
configuration-averaged 
(DCA).  Previously  known  as 
electron  configuration 


]  LANL 
I  DCA 
■LANL 
■  DCA 
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Argon  CR  Test  Cases 


Isothermal  Heating 


•  Comparison  between  DCA,  QSS, 
and  single  Boltzmann  grouping 

•  Single  Boltzmann  grouping 
applies  grouping  to  whole  ionic 
distribution 

•  "Worst  case"  scenario  for 
accuracy 

•  Most  comparable  formulation 
to  QSS 


Irradiated  System 


•  Implemented  radiation  source 
terms  for  the  DCA  case  with 
Planckian  field  (i.e.  black  body 
radiation) 

•  Required  addition  of  electron 
energy  equation 

•  Uses  tabulated  cross  sections  to 
calculate  photo-ionization  rates 
(more  flexible  than  analytical) 
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Test  Case:  Isothermal  Heating 


Conditions 

•  Heavy  ions: 

—  1025  rrr3 

•  Initial  electron  density 

— 1010  rrr3 

•  Fixed  electron 
temperature 

-50  eV 


CO 


a 

TJ 

o 

+3 
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0) 
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Atomic  <Z> 


Test  Case:  Isothermal  Heating 


Conditions 

•  Heavy  ions: 

-  1025  nrr3 

•  Initial  electron  density 

-  1010  nrr3 

•  Fixed  electron 
temperature 

-  50  eV 


Energy  Level  [eV] 


23  October  2017 
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Test  Case:  Isothermal  Heating 


Conditions 

1035 

1030 
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1020 

1015 

icr12  to-11  io-10  icr9  kt8 

Time  [secs] 
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Heavy  ions: 
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Test  Case:  Irradiated  System 


Conditions 

•  Heavy  ions: 

-  1025  nrr3 

•  Initial  electron 
density 

-  1010  nrr3 

•  Fixed  radiation 
temperature 

-  50  eV 
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Atomic  <Z> 


Test  Case:  Irradiated  System 


Conditions 

•  Heavy  ions: 

50 

-  1025  nr 3 

•  Initial  electron  _  40 
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Test  Case:  Irradiated  System 


Conditions 

•  Heavy  ions: 

-  1025  nrr3 

•  Initial  electron 
density 

-  1010  nrr3 

•  Fixed  radiation 
temperature 

-  50  eV 


Energy  [eV] 
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Non-Maxwell  an  Inelastic  Collision  Modeling:  Overview 


*  Model  inelastic  and  elastic  collisions  amongst  electrons,  hydrogen  atoms,  and  ions  while  allowing  non-Maxwellian 
electron  energy  distribution  functions.  (Excitation,  De-excitation,  Ionization,  Recombination,  the  Fokker-Planck 
Equation,  Bremsstrahlung,  along  with  radiative  transitions  and  many  other  processes  in  the  future). 

9  Previous  work  on  this  task  either  assumed  a  Maxwellian  electron  distribution  function  OR  involved  a  very  simple 
collision  model  (BGK). 

9  The  task  objective  is  not  the  further  development  of  these  physical  process,  or  even  analytical  approximations  of 
them,  but  rather,  an  attempt  to  incorporate  as  much  of  the  important  physical  process  as  possible  with  as  little 
computational  cost  as  possible. 

9  To  that  end,  we  are  in  the  early  processes  of  characterizing  all  the  different  types  of  relevant  CR  processes  that 
occur,  along  with  the  assumptions  on  shape  of  the  distribution  functions,  and  determining  the  circumstances 
under  which  higher  accuracy/higher  computational  cost  models  are  necessary  to  fully  implement. 
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Non-Maxwellian  Inelastic  Collision  Modelling 


•  Complex  CR  Formulas: 

—  Typical  differential  equations  can  be  written  on  a  single  line  where  as  these  take  pages  to 
write  out. 

•  Boundary  conditions: 

—  Vlasov  boundary  conditions  are  a  challenge.  Additionally,  each  inelastic  collision  process 
has  its  own  type  of  boundary  condition  in  which  a  region  can  be  drawn  out  where  there  is 
sufficient  energy  for  that  process  to  occur. 

•  Non-local  interactions: 

—  Even  for  explicit  methods  with  a  small  time  step,  updating  the  solution  at  a  single  location 
requires  calculations  with  a  dependence  on  every  point  being  modeled. 

•  Incorporation  of  other  non-trivia  I  mathematical  algorithms: 

—  Pre-computation  of  quantities  required  for  propagation  of  the  solution. 
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Non-Maxwellian  Inelastic  Collisions:  Test  Case  Setup 


*  Atomic  Hydrogen 

*  Electron  impact  excitation/de¬ 
excitation 

*  40  excited  states 

*  Delta  function  for  H,  H+ 


*  Photo-ionization  (pre-process) 

-  41  nm  laser 
-30  eV 

*  Initial  Conditions 

—  Atmospheric  like  conditions 
—  H  =  2.68xl025nrr3 

-  H+,  e  =  4.06xl021  rrr3 

-  Te  =  10.9eV 
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Normalized  Electron  Density 


Non-Maxwellian  Inelastic  Collision  Modeling:  Laser  Example 

(preliminary) 
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Boltzmann  Ratio  of  two  highest  energy  species  (for  the 
same  simulation  shown  on  the  left) 
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Initial  photo-ionization  of  hydrogen.  Discrete  changes  in  electron  velocities  from 
excitation/de-excitation  and  smoothing  from  the  Fokker-Planck  Equation 
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Non-Maxwell  an  Inelastic  Collision  Modeling  :  Limitations  of 

Standard  Numerical  Methods 


•  Multiscale  Hybridization: 

—  Standard  explicit  schemes  may  not  perform  adequately  when  modeling  highly  non¬ 
polynomial  solutions  and/or  when  the  number  of  points  employed  to  represent  the 
solution  is  huge.  An  approach  to  improving  performance  is  to  fall  back  on  a  known 
low  computational  cost  coarse  solution  approximation  rather  than  a  polynomial. 

•  Sobolev  Regularization  with  Coordinate  Transforms: 

—  Implicit  schemes  may  be  less  than  desirable  when  they  require  solving  complicated 
non-linear  equations  where  the  repeated  linear  solves  to  be  performed  involve 
dense  matrices.  Previous  work  has  shown  this  type  of  regularization  to  be  an 
effective  alternative  to  fully  implicit  methods  for  complex  problems. 

•  Lagrangian-TV  Regularization : 

—  Another  type  of  regularization  that  may  better  allow  for  shocks  and  other  specific 
solution  features. 
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Regularization:  Motivation 


•  Faster  Linear  Solves 

—  The  matrices  that  arise  from  implementing  the  implicit  schemes  for  the  differential  equations  that  we  wish  to  solve  are  dense 
and  ill  conditioned.  The  matrices  arising  from  the  regularization  are  always  sparse  and  well  conditioned. 

•  Faster  Development 

—  The  differential  equations  are  complex  and  an  implicit  method  needs  to  incorporate  all  aspects  of  /  (boundary  conditions, 
convolutions,  precomputing  global  variables,  etc.)  which  can  be  a  time  consuming  process. 

•  Vlasov  Merging 

—  Lagrangian  coordinates  with  frequent  re-gridding/re-initialization  (semi-Lagrangian)  are  a  standard  approach  for  modelling  the 
Vlasov  equation.  Looking  ahead  to  a  Vlasov  simulation  with  inelastic  collisions  and  other  non-linear  processes,  it  may  be 
desirable  to  have  a  single  Lagrangian  framework  to  model  all  these  processes. 

•  Lagrangian-7V Regularization  on  u(X)  has  a  'light  touch ' 

—  In  the  extreme  case  where  the  time  step  is  small  enough  that  no  regularization  is  needed,  the  solution  will  not  change  at  all. 

•  Two  Step  Regularization 

—  Even  in  cases  where  a  smooth  solution  is  sought,  applying  the  TV-Regularization  first  ensures  even  a  small  amount  of  L2 
smoothing  will  bring  the  solution  into  a  Sobolev  space. 
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Regularization:  General  Background 


The  standard  approach  to  achieving  stability  of  a  numerical  scheme  is  to  write  it  in  such  a  way  that  the  propagation  in 
time  is  implicitly  defined.  A  pen  and  paper  analysis  of  the  scheme  (like  Von  Neumann)  is  used  to  establish  stability, 
which  typically  means  that  the  function  will  remain  in  a  Sobolev  space. 

•  Explicit:  zn+1  =  zn  +  dtf(zn )  Implicit:  zn+1  =  zn  +  dtf(zn+1 ) 


An  alternative,  somewhat  more  recent  and  less  common  approach  is  to  keep  a  solution  in  a  Sobolev  space  by  simply 
directly  incorporating  that  restriction  into  the  numerical  method: 


yn+1=argminj(y-dtf(zn)ydX  +  j(VyydX 


Zn+1  =  z„  +  yn+i  =  Zn  +  0  —  A)  1(dtf  (zn)) 


While  extremely  simple,  this  approach  has  been  shown  effective  for  even  very  complex  simulations. 

An  alternative  functional  space  to  consider,  most  common  in  data  processing  but  increasingly  spreading  into  physics 
applications,  is  the  following  (with  one  L2  term  and  one  L1  term): 

•  yn+i  =  argminj (y  —  dtf(zn ))2  dx  +  f\Vy\dx 


And  more  recently,  fast  approaches  to  solving  problems  like  this  that  do  not  rely  on  classical  approaches  (i.e.  gradient 
descent  and  the  Euler-Lagrange  equation)  have  been  developed. 

Regularization  can  be  applied  directly  to  a  distribution  cp(x )  (Eulerian)  or  to  u(X)  =  x  (Lagrangian). 
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Lagrangian-TV  Regularization:  Burger's  Equation  Example 


•  The  Lagrangian-TV  regularization  can 
produce  an  exact  shock  solution  with  no 
diffusion  and  no  dispersion  in  a 
computationally  efficient  and 
straightforward  manner. 

•  1st  and  3rd  Order  solutions  provide 
accurate  shock  locations  but  have 
diffusion  or  dispersion  errors 

•  Lagrangian-Sobolev  does  not  get  the 
correct  shock  location  and  has  moderate 
diffusion  errors 

•  Euler-TV  and  Euler-Sobolev  does  not 
obtain  the  correct  shock  location  and  has 
significant  errors  in  its  profile. 


Inviscid  Burger's  Equation,  (pt  =  (p(px,  solved  with  various  schemes 
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Conclusion 


*  Maxwellian  CR 

—  Added  QSS  and  single  temperature  Boltzmann  grouping  and  compared  with  DCA 
results 

—  Included  radiation  source  term  with  Planckian  field  which  is  a  stepping  stone  to 
radiation  transport 

*  Non-Maxwellian  CR 

—  Demonstrated  ability  to  capture  non-Maxwellian  EEDF  effects  on  Zero-D  LPI  like 
simulation 

—  Implemented  the  Fokker-Planck  collision  operator  along  with 
excitation/dexcitation 

—  Ionization  and  recombination  processes  are  near  completion 


23  October  2017 
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Future  Work 


*  Maxwellian  CR 

—  Affect  of  grouping  on  generated  spectra 

—  Automatic  level-grouping  (looking  for  collaborators!) 

—  Coupling  with  fluid  model  for  ID  simulations  (high  Mach  Argon  shock  ID  LPI 
breakdown) 

—  Use  cFAC  to  generate  atomic  data  for  low-  to  mid-Z  elements  (up  to  Xenon) 

—  Compare  with  other  CR  codes,  e.g.  FLYCHK 

*  Non-Maxwellian  CR 

—  Add  Bremsstrahlung,  along  with  radiative  transitions 

—  Move  to  2V  Fokker-Planck  solver 

—  Further  implement  regularization  techniques  to  collisional  modelling  and 
compare  against  standard  explicit  and  implicit  schemes 
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Questions? 
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